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FROM BOSONIC GRAND-CANONICAL ENSEMBLES TO NONLINEAR 

GIBBS MEASURES 

NICOLAS ROUGERIE 


Abstract. In a recent paper, in collaboration with Mathieu Lewin and Phan Thanh 
Nam, we showed that nonlinear Gibbs measures based on Gross-Pitaevskii like functionals 
could be derived from many-body quantum mechanics, in a mean-field limit. This text 
summarizes these findings. It focuses on the simplest, but most physically relevant, case 
we could treat so far, namely that of the defocusing cubic NLS functional on a ID interval. 
The measure obtained in the limit, which lives over H 1 ^ 2 ~ e , has been previously shown to 
be invariant under the NLS flow by Bourgain. 
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1. From many-body quantum mechanics to nonlinear Schrodinger theory 

The nonlinear Schrodinger equation is a basic model of quantum optics and condensed 
matter physics, whose analysis has motivated a huge mathematical literature. The equa¬ 
tion itself has been originally introduced as a phenomenological model, and its derivation 
from the fundamental axioms of quantum mechanics is a major issue in theoretical and 
mathematical physics. 

The advent of Bose-Einstein condensates of ultracold dilute alkali vapors in the mid- 
90’s has given a renewed impetus to the field. The NLS equation is nowadays thoroughly 
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studied as the basic effective model for the description of such objects. The reason why a 
one-body model such as NLS theory can be used to describe a many-body system is that, 
in this context, all particles reside in the same quantum state (this is what is meant by 
Bose-Einstein condensation). 

An important line of research has been the rigorous derivation of the equation, in a certain 
limit, from the many-body Schrodinger equation for bosonic particles. This section recalls 
some background on this topic. The author apologizes for not attempting to do justice to 
the relevant literature in any exhaustive way: the reader is refered to 13 |22J m, HU M S3] 
for a more complete introduction to the subject. 

1.1. Derivation of ground states. The basic object one starts from is a many-body 
Schrodinger Hamiltonian, say of the form 

N 

h n = y (~a j + v ( Xj )) + y w N (xi-xj). (i.i) 

j =1 1 ^i<j^N 

Here V is an external potential felt by the particles and wjy a pair-interaction potential, that 
is most often taken radial. This operator is considered as acting on L 2 (R dN ) ~ (^)( v L 2 (R d ), 
the space of square integrable, symmetric A r -body bosonic wave functions, that is of those 
% G L 2 (R dN ) ~ L 2 (R d ) satisfying 

^N 1 1 • • • i %N ) ^N (^cr(l) ? * * • ? ■^cr(A^)) 

for all permutation cr of the coordinates x\,... ,xn £ R d of the particles under consideration. 
This symmetry requirement means that we consider bosonic particles, which, as the name 
indicates, is mandatory to obtain a Bose-Einstein condensate. On a physical level this 
implies that we consider the atoms of the gas as the basic particles, ignoring that they are 
themselves made of more fundamental objects. This is however as fundamental as one can 
get for a reasonble description. 

By ground state we mean the state of lowest possible energy, i.e. the eigenfunction of the 
operator achieving the lowest eigenvalue, or, equivalently, solving the minimization problem 
for the quadratic form 

£n[^n] = ( 1 . 2 ) 

under the constraint that Tjv € L 2 (R dN ) be .^-normalized. This gives the most stable 
equilibrium configuration at zero temperature. 

An important question is the following: is it true that, at least in some well-defined scaling 
limit, the ground state(s) of the above Hamiltonian resembles a Bose-Einstein condensate? 
The latter is a IV-body wave function of the form 

N 

ty N (xi,...,x N ) = JJrt(xj), (1.3) 

3 =1 

corresponding to putting all the particles in the same quantum state u £ L 2 (W l ). In classical 
mechanics this would correspond to assuming independent identically distributed particles, 
i.e. molecular chaos, Boltzmann’s Stosszahlansatz. 
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If such an ansatz is valid, one naturally obtains a NLS-like energy functional that u 
should minimize. Indeed, inserting in (11.211 we find 

= N ( [ l Vu | 2 + + —^—7 ,—— [[ \u(x)\ 2 w n (x - y)\u{y)\ 2 dxdy. 

\J K d / 1 JJ R d xR d 

One usually expects the ansatz (11.311 to be valid for large particle number N. In order that 
the one- and two-body terms in the above functional weigh the same in this limit, so that 
one might hope for a well-defined limit, it is natural to assume that wn be of order IV — 1. 
At least two possibilities to achieve this should be considered: 


• For a fixed potential w, pick 


w n 


A 


IV — 1 


w 


for some A € M kept fixed (or converging to a constant) when N —>• oo. This situation is 
often called the mean-field limit, it physically corresponds to having many weak collisions 
between the particles. The limit functional one obtains is then 


S[u] = [ |Vu| 2 + V\u\ 

Jm d 


+ 2 


\u(x)\ 2 vj(x — y)\u(y)\ 2 dxdy , 


(1.4) 


i.e. a nonlinear functional with non-local nonlinearity, often called Hartree’s functional. 


• For a fixed constant /3 > 0 and a fixed potential w, pick 

w N = v A ^ N d ^w(N^.) 

with A as above. This can be called a NLS limit. It allows to take the range of the potential 
as a parameter, and to describe a situation where, because of diluteness, collisions can be 
rare but quite strong. This happens when /3 > 1 /d. Since 

N dfi w{N< 3 .) (J d *>) 6 o 


as measures, it is natural to expect a true NLS functional with local nonlinearity 


£[u]= f |Vu| 2 + V\u\ 2 + ^ f f w\ f \u(x)\ 4 dx (1.5) 

JR d 2 \J^d J 


in the limit. 


We refer to the literature, in particular the previously mentioned reviews [331141] and 
also m mm na [39] for results in this direction. Under fairly general assumptions, one 
may obtain convergence of the energy per particle: 

E(N) 

——- > e when A —>• oo 

N 

where E(N) is the bottom of the spectrum of and e the infimum (under a unit L 2 - 
mass constraint) of the relevant limit functional. This serves as a justification for the use 
of NLS functionals to describe the ground states of large bosonic systems. We warn the 
reader that for f3 relatively large in the NLS limit, important new aspects emerge and an 
ansatz (11.31) is no longer sufficient, see [32l 136] as well as [33] and references therein. The 
most prominent case where these important subtleties occur is when d = 3 and [3 = 1, the 
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so-called Gross-Pitaevskii limit. Then, short-range correlations in the ground state have 
the effect of replacing f R3 w by 4-7ra in (11.51) . where a is the scattering length of w. 

One can also prove results on (sequences of) ?V-body minimizers, showing that in a certain 
specific sense, they are close to factorized states of the form (11.31) . Let us stay vague on this 
point for the moment and simply warn the reader that in general, the wave function Tv 
itself is not close to a condensate, only its reduced density matrices are. 


1.2. Derivation of evolution equations. Once one knows that ground states of large 
bosonic systems show Bose-Einstein condensation, the next natural question is wether this is 
preserved by the many-body Schrodinger dynamics. Indeed, a typical experimental situation 
is as follows: create a Bose-Einstein condensate by cooling a Bose gas to extremely low 
temperatures, so that one may assume that the ground state is reached. The IV-particles 
are then in a stationary state of the Hamiltonian Hn which is (approximately) a factorized 
state u® N . To probe the dynamics one may then perturb the system, most often by changing 
the one-body potential V in (11.11) and observe the subsequent evolution. 

This asks the question of the dynamics of (approximately) factorized states under the 
quantum evolution associated to Hamiltonians of the form (HU- Consider the Cauchy 
problem 

[idt^N = H n ^n ,, 

i'Mcn^r (L6) 

where we stay voluntarily vague on the precise meaning of the symbol ~. Usually this 
means convergence of reduced density matrices, the typical result one obtains for ground 
states. Is it true that the solution fjv(l) stays approximately factorized 

N 

^ N (t,X 1 , . . .,X N ) « ( l7 ) 

3 = 1 

after some time evolution? In limits similar to those previously mentioned, this has been 
proved by a variety of authors. Inserting the ansatz in (11.61) one sees that the relevant u(t) 
solves a nonlinear equation: 


• Hartree’s equation 


idtu = —A u + Vu + X(w * |it| 2 )u 
u(0) = Uq 


in the case of a mean-field limit. 


( 1 . 8 ) 


• a NLS equation 


idtu = —A u + Vu + (f Rd w) A|u| 2 u 

u(0) = Uq 


(1.9) 


in the case of a NLS or Gross-Pitaevskii limit. In the latter case, / Rd w should be replaced 
by 47 tx (scattering length of w). 


Many results in this direction have been proved over the years, see the reviews [3 EOS! 
and for example the (by far non exhaustive) following list of original references [I] (23# 2D, 

na m nasoi Eg • 
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1.3. Positive temperature stationary states. The results quickly summarized in the 
previous subsections already put the NLS theory on a firm rigorous ground as far as the 
description of Bose-Einstein condensates at low temperature is concerned. But it would 
be highly desirable to know something about not only ground states as in Subsection 11.11 
but also positive temperature states. Indeed, in experiments the temperature is always 
(extremely) low but finite. The Bose-Einstein condensation phenomenon in fact refers to 
the existence of a critical temperature below which the positive-temperature equilibrium 
states of a Bose gas show condensation in the sense that all particles reside in the ground 
state of NLS functionals as previously defined. A satisfactory, mathematically rigorous, 
estimate of the critical temperature is still missing. 

A closely related question is that of the description of stationary, positive temperature, 
states in the absence of Bose-Einstein condensation. Can the NLS description still be 
relevant, in a well-defined limit, for such objects? The (canonical) equilibrium state for the 
Hamiltonian at temperature T is the Gibbs state 

= £—exp (1-10) 

a trace-class operator on L^(M. dN ), where the partition function Zv,t fixes the trace to be 1. 
Can such states, or close variants, converge to NLS-related objects in some limit? What 
would the natural limiting objects be? These questions are the subject of the paper m 
that we summarize in this text. Here we shall first try to answer the second question. Since 
many-body Gibbs states are stationary under the many-body Schrodinger flow, one should 
expect that the natural limit objects are invariant under the NLS flow. 


2. Nonlinear Gibbs measures and the NLS flow 

Nonlinear Gibbs measures have recently become a useful tool to construct solutions to 
time-dependent nonlinear Schrodinger equations with rough initial data, see for instance [261 
mm EDi nans™ in]- These are probability measures which are formally defined by 

dn(u) = ^e~ s ^ du, (2.1) 

Zj 

where 8[u] is a NLS-like functional as discussed above. The rigorous construction of these 
measures is an important part of the Euclidean approach to Constructive Quantum Field 
Theory (CQFT) [161 21, 33 03 13 SI]- The recent progress on the PDE side has mainly 
been to define the NLS flow almost surely on the support of the measures, which typically 
live over rather rough functional spaces, and to prove the invariance of the measures along 
the associated NLS flow. 

These objects have an obvious formal resemblance with Gibbs states such as (11.1011 and 
it is the main goal of the paper m to explore this connection in more details. For the 
moment we can treat only the defocusing ID case based on taking d = 1 and A, w ^ 0 
in (HD or m- This is due to well-known difficulties with the definition of the measures 
in other cases. For this reason the rest of the discussion will be mainly limited to the case 
d = 1, although we have partial results in a more general and abstract setting, for which 
we refer to the original paper ED- 
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2.1. A reminder on classical dynamics. Consider the evolution of a single classical 
particle. The dynamics is given by a Hamilton function H(x,p ) : M d xl rf ^Ion the phase 
space of all possible positions and velocities. The Hamilton equation for the dynamics of 
the particle is 

(x = V p H(x,p) 

= -V x H(x,p). 

Define the associated Gibbs measures at temperature T 

dp(x,p ) = — exp (— T~ 1 H(x,p)) dxdp (2-2) 

Zj 

where Z is a normalization factor ensuring 


dp(x,p) = 1. 


These measures are of course invariant under the flow, because the latter conserves the 
energy H{x,p ) and leaves the Lebesgue measure invariant (Liouville’s theorem). One should 
note that the Gibbs measures are not any invariant measure, in that they have a physically 
very natural variational characterization. Namely, they minimize the free-energy functional 


H(x,p)dp(x,p) + T 


p log p dxdp 


amongst probability measures, which means that they achieve a balance between minimizing 
the energy and maximizing the entropy. 

Again by formal analogy, between m and (12.2|) this time, the measures m are good 
candidates to play a similar role in the NLS theory where the dynamics is (at least formally) 
also Hamiltonian. 


2.2. Invariant measures for the NLS equation. Let us now recall the rigorous defini¬ 
tion of the Gibbs measure m , based on the functionals El or El- The way to proceed 
is by first defining the free, non-interacting, measure (case A = 0), formally 

1 


po{du) = — exp — / |Vitp du 
\ Jn J 

where Q C K rf . Then one defines the interacting measure as being absolutely continuous 
with respect to pq: For a given function w or even for w = 5q, set 


d(du) = -Jr exp 

Zj r 


flxfl 


\u(x)\ 2 w(x - y)\u{y)\ 2 p 0 (du ) 


(2.3) 


where the number Z r (relative partition function) ensures that this is a probability measure. 
In [27] we have considered only cases where this construction makes sense without any need 
for renormalization. 
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The free measure. The free measure is a gaussian measure over an infinite dimensional 
Hilbert space. It is defined by the following classical construction [8l 06] . Take h ^ 0 a 
self-adjoint operator on some Hilbert space S) with compact resolvent. For the applications 
in this note we shall take 

d 2 

+ |x|“ + m (2.4) 


h = - 


dx 2 


on 


or 


d 2 

h = ——r -I- m 
dx z 


(2.5) 


on [—1,1] C M with either, Dirichlet, Neumann or periodic boundary conditions. The L 2 
space on which these operators act will always be denoted 1} = Lr(M). We require that 
m € R is such that these operators are positive definite. In the latter two cases one must 
pick m > 0 to ensure this. We shall restrict to a > 2 for reasons that we will explain later, 
and include the case (1231) as a formal limit a —> oo of m- 
We write the spectral decomp ositiorQ 

OO 

h = J2 x j\u j ){u j \ 

3 = 1 

and define the asociated scale of Sobolev-like spaces 


f) s := D(h s / 2 ) = <j u = ^ i 

m 


\u\\l* 


m 


( 2 . 6 ) 


One can then define a finite dimensional measure on span(ui,..., uk) by setting 


K 


A 


ddo(u) := (^) -^exp (-Xj\(u, Uj)\ 2 ) d(u,Uj 


j= 1 


where d{u,Uj) = dajdbj and aj,bj are the real and imaginary part of the scalar product. 
A general tightness argument for this sequence of measures allows to see that they are all 
cylindrical projections of a common measure living on some Sj s , where generically s < 1: 


Lemma 2.1 (Free Gibbs measures). 

Assume that there exists p ^ 0 such that 

M = 


3=1 3 


< OO. 


(2.7) 


Then there exists a unique measure po over the space $) l ~ p such that, for all K > 0, the above 
finite dimensional measure po t K is the cylindrical projection of po on span(ui,..., uk)- 
Moreover 


7( W : = / \ u ® k ){u® k \ dp 0 (u) = k\ {h^f 

Jfp-p 


( 2 . 8 ) 


where this is seen as an operator acting on 
Some comments: 


^jv)(t>| denotes the orthogonal projector onto v € Sj. 
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• In the ID cases with anharmonic potentials that are our main example, one can pick 
p = 1, so that the measure lives on the original Hilbert space L 2 (M d ). In general 
this is not the case, and this is the main reason why we have so far been able to 
derive the NLS Gibbs measures only in ID. 

• One can check that in the case of a finite interval (a = oo formally), (12.71) holds 

for any p > 1/2. The measure thus lives over 1,1]), the usual 

L 2 -based Sobolev space, for any e. More generally, for a Laplacian on a bounded 
set in d dimensions, the measure lives over H 1 ~ d / 2 - £ for all e > 0. We thus see that 
already in 2D, one has to consider measures over negative Sobolev spaces. 

• One might be surprised that the measure lives over a space where the energy asso¬ 
ciated with the Hamiltonian h is infinite. This is however a basic fact in the theory 
of gaussian measures. 

• Formula (12.81) indicates that although the measure itself may live over rough spaces, 
certain averages related to it are in fact nice compact operators over the original 
Hilbert space. In the cases of interest to this note, the right-hand side of (12.81) will 
actually always be trace-class. 


The interacting measure. We can now use the previous lemma to define the interacting 
measure. We have to check that 


-FnlM := 



y)\u(y)\ 2 dxdy 


makes sense /zo-almost surely and satisfies some integrability properties w.r.t. [ 1 q. We only 
need that the relative partition function Z r makes sense and is strictly positive, but in the 
simple case under consideration we will have much stronger properties. 


Lemma 2.2 (Interacting Gibbs measure). 

Consider the case where h is either (E3D or m, let p be the number in m and po the 
free measure defined in Lemma FOl Assume that w = w\ + W 2 where w\ is a positive Radon 
measure and W 2 a positive bounded function. 

The function u i->- FnlM is in L 1 dpfi). In particular 

p{du) = exp (—TnlM) Po(du) 

/LJ y 

makes sense as a probability measure over . 


The proof consists of essentially two remarks: 

• Since we work with a positive w, it is sufficient to check integrability properties for 
Fnl[«] itself to obtain Z r > 0. This makes computations easier. 

• That FnlM is integrable follows easily from (12.81) . Identifying w with the multipli¬ 
cation operator by w(x — y) on the two-body Hilbert space, we have 

J F NL [u\dp 0 {u) = \ J Tr [ w |tt® 2 )(«® 2 |] dp 0 (u) 

= j Tr [w h~ x ® h~ l ] , 

where the trace is taken over the symmetric subspace. The expression in the right- 
hand side can be shown to be finite by direct estimates on the Green function /i _1 . 
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When w is just a bounded function one can directly use the fact that h 1 is trace- 
class. 

Invariance under the NLS flow. For the particular ID case (12.51) . Bourgain [9] proved that 
the NLS flow is well-defined /i-almost surely and that the measure is invariant under the 
flow0. He also considered higher order non-linearities, as well as the focusing case. Still in 
ID, the harmonic oscillator case (a = 2 in m, not covered by our approach) has been 
dealt with by Burq-Thomann-Tzvetkov [12] . Here also one gets existence of the flow p-a.s. 
and invariance of the measure. 

These measures thus play a special role in the theory of the NLS flow: they give a 
whole set of low regularity initial data for which the flow makes sense, and they are a 
natural invariant for the dynamics. Since the NLS equation arises as the limit of many- 
body quantum mechanics, it is very natural to wonder wether these objects can be related 
to equilibrium states of some many-body Schrodinger Hamiltonian. We explain that this is 
indeed the case in the next section. 

3. Derivation of the NLS-Gibbs measure on a unit interval 

In this section we turn to describing the new results proved in m ■ For simplicity we 
shall not present the full abstract setting but only those results that are relevant for the 
models based on (12.41) and (12.51) . for which (12.71) holds with p = 1 and Lemma 12.21 holds 
trucjE It is in fact possible to relax these assumptions and obtain a measure which lives 
over negative regularity spaces, whose density matrices (12.81) are not trace-class but only 
belong to a higher Schatten space. In this case the interaction operator w must be very 
regular and cannot be of the form w(x — y). One should think of it as a regularization of a 
physical interaction. 

The many-body object that gives rise to the nonlinear Gibbs measure is in fact the grand- 
canonical ensemble for bosonic particles. The limit we consider is a large temperature/mean- 
field limit. We shall not try to discuss the physical relevance of this limit here. The main 
goal is to connect nonlinear Gibbs measures to many-body equilibrium states. In view of 
the results recalled in Section Q] it is then very natural to expect the former to be invariant 
under the NLS flow. 

3.1. Grand-canonical ensemble. The grand-canonical Gibbs states live over a Hilbert 
space where the particle number is not fixed, but seen as a random variable. Namely, 
we work in the bosonic Fock space made of all possible superpositions of states with any 
number of particles: 

T = C ® Sj © S) 2 © ... © Sy N © ... 

T = C © L 2 (D) ffi L 2 (H 2 ) © ... © L 2 S (Q N ) © ... 

where !lcR and Sj ) = L 2 (H). Everywhere in the sequel we use the notation 

k 

:=03 

S 


2 Stricto sensu he worked under periodic boundary conditions only. 
^This is refered to as “the trace-class case” in m■ 
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for the fc-particles bosonic Hilbert space. In this context the particle number is an operator, 
acting on each n-particle sector separately, and equal to nl on the n-particle sector: 

oo 

Af : = n 11 

n =0 

The mixed states over this Hilbert space are as usual the self-adjoint, positive, trace-class 
operators on T with trace 1. The pure states are the vectors of J-, identified with the 
corresponding orthogonal projectors. For any operator A on T we call its expectation value 
in the state T the quantity 

(A) r ■= Trjr [AT ], 

whenever this has a meaning. In particular, the expected or average particle number is 

Trjr [NT]. 

We consider the second quantized Hamiltonian acting on J 7 , which is the natural extension 
of (ED: 

e ° = © (j = © (+ \ x i\ a+m j 

00 / \ 00 / 

w = © I] w *3 = ® I] w ( x i - x i) 

n= 2 yl J n= 2 yl 

00 

H a = Ho + AW = ®^„ )A . 

n= 1 

Here h is taken as in (12.41) or (12.51) (which we include as the formal a = oo case). The 
interaction potential w is chosen as in Lemma 12.21 in particular w ^ 0, and we introduced 
a coupling constant A ^ 0. Note that this choice of w is infinitesimaly form-bounded with 
respect to h , even if w = So, because of the ID Sobolev embedding. The Hamiltonian thus 
makes sense as a self-adjoint operator. 

The Gibbs state at temperature T for this Hamiltonian is by definition the minimizer of 
the free-energy functional 

T A [T] = Trjr [H a T] + T Tr^ [T log T] (3.1) 

defined for mixed states over T. The first term is the energy, the second the von Neumann 
entropy. Minimizing over all states we clearly obtain 

rA ’ T= ziM exp ( _ ^ Hi ) (3 ' 2) 

where the partition function Z\(T) fixes the trace equal to 1 and satisfies 

F X (T) = -Tlog Z X (T) 
where F\(T) is the minimum free energy. 


(3.3) 
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It is useful to introduce the notion of fc-particle reduced density matrix for states oved 
T. For a “diagonal” state of the form 


r — G n , G n : S) n t—>■ ij n , 


n =0 


this is the operator on the fc-particles space Sj k dehned b}0 


Tr f 


rCzk 


A k r^] = ]T (J Tl '^ W 10 {n ~ k)G 


(3.4) 


for every bounded operator A k on S) k . These are helpful for rewriting expressions involving 
only fc-particle operators. In particular the number operator involves only one-body terms, 
so the expected particle number can be computed as 

pi 1 ) 


Ti> [NT] = Ti s - 


-r r J 

The energy involves only two-body terms, and we have: 


TY^pHUr] = Tr^i [hr (1) ] + A Tr ^2 \w T (2) 


(3.5) 


(3.6) 


The limit we shall consider is 

T —»■ oo, AT->1. (3.7) 

The rationale here is that to see a non-trivial effect of the temperature in a bosonic system, 
T should be rather large, because of the Bose-Einstein condensation phenomenon we alluded 
to before. Then one can easily get convinced that the expected particle number behaves as 
T in this limit. For example, consider the A = 0 case. The one-particle density matrix can 
be computed explicitly: 


■p(i) _ 

1 0 ,T ~ 


1 


e h/ T _ x' 

After dividing by T, this clearly converges to the k = 1 case of (12.811 . Since we assume the 
latter object to have a finite trace, in view of (13.51) . it is clear that the expected particle 
number of the free Gibbs state grows like T. Taking the coupling constant A ~ 1/T is thus 
the equivalent, in the grand-canonical setting, of the canonical mean-held limit discussed 
in Section Q] where the coupling constant scaled as the inverse of the particle number. Note 
that in case p > 1 in dm this argument is no longer valid, and the particle number will 
grow much faster than T. This makes this case much more difficult to deal with, which is 
why it is not included in this note and we refer to the original paper. 


3.2. Main theorem. Let us now state the main result of this note, proved in [27]. It 
concerns the limit (13,7p of the previously dehned grand-canonical ensemble. One should not 
expect a convergence of the Gibbs states themselves since their expected particle number 
diverges. Neither should one hope for an estimate in some norm, for well-known reasons 
that we shall not elaborate on here. The correct objects to look at are the £;-body density 
matrices with k much smaller than the particle number (in particular, k hxed will do). 
The results on ground states and dynamics we have mentioned in Section Q] are all of the 
form of estimates on such objects, and so is the following theorem on positive temperature 
equilibrium states: 

4 The partial trace is taken over the symmetric space only. 
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Theorem 3.1 (Derivation of nonlinear Gibbs measures). 

Let h be given by (12.41) or (12.51) . and w satisfy the assumptions of Lemma \2.2l Let p be the 
Gibbs measure defined therein. For every fixed k 1 we have 

J^\u« k )(u« k \d»(u) (3.8) 

strongly in the trace-class norm, in the limit CUD. 


Some comments: 

• When L' 2 (Ll) is replaced by a finite dimensional Hilbert space, canonical versions of 
this theorem have been proved before, see [23l [25] and m Appendix B]. 

• The collection of all the “moments” of p for k = 1,2... appearing in the right-hand 
side of (|3.8I) characterizes p uniquely. This measure is thus the natural, unique, 
limit object for the grand-canonical ensemble. 

• Roughly speaking, (13.81) is a way of making rigorous the approximation 

rA,T~ J \ti(VTu))(aVfu)\dp(u) (3.9) 

where £(u) is the coherent state built on v : 


€(v) 


= e 


IMI 2 /2 





v ®k 

7H‘ 


(3.10) 


3.3. The variational formulation. The proof of Theorem [XII is variational. The conver¬ 
gence of density matrices is obtained as a corollary of estimates on the minimum free energy, 
or equivalently, the partition function. An important difficulty is however that Z\{T ) is dif¬ 
ficult to control directly, in particular because even at A = 0 this quantity diverges very fast 
when T —>• oo. Our strategy is therefore to instead deal with the difference between the free 
energies at A = 0 and A > 0, that is, to count everything relatively to the free Gibbs state. 
This is possible because, as one easily realizes, also minimizes the relative free-energy 
functional 


T iel [T] = A Tr T [W T] + TTr T [T (logT - log T 0 ,t)] 


= A Ti’^2 



+ T TV j- [r (log r- log To,r)] 


(3.11) 


amongst all mixed states. The second term in the above is the (quantum) relative entropy 
of r with respect to ro,T> which is always positive. The infimum of the above functional 
equals 

f x (t)-Fq(t) = riog^j^y. 

Studying the minimization of (13.111) thus allows for a direct access to the ratio of the inter¬ 
acting and non-interacting partition functions, which will converge to some limit although 
the partition functions themselves diverge. We obtain Theorem 13.11 as a corollary of the 
estimate 
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Theorem 3.2 (Asymptotics of the free-energy). 

Let Z r be the relative partition function for the Gibbs measure fi, defined in Lemma \2.2l 

Z r = J exp (-F N lM) df-io(u) 

with no the free Gibbs measure. In the limit (EH) we have 


F X (T) - Fp(T) 
T 


—y — log Z r . 


In order to complete the variational framework of the proof, one characterizes — log Z r 
as the inhmum of some variational problem of which p is the solution. Indeed, one easily 
sees that [i minimizes 

do 
dp o 

= J I [w |«® 2 )(u® 2 |] <M») + J±- log (JFj i„ a (3.12) 

amongst all /zo-absolutely continuous measures v on the one-body Hilbert space S). The 
second term is the classical relative entropy of u relative to ho- 

The semi-classical nature of Theorems im and is now apparent. We have to relate 
the minimization problems for the quantum and classical free energy functionals (13.111) 
and (13.121) . This means that we need to understand how to describe quantum objects, 
operators on a Hilbert space, in terms of classical objects, probability measures. This 
requires specific tools, and we shall review some of these in the next section. 

Since we work on a variational problem, we can provide separately an upper and a lower 
bound on the free energy. The former essentially follows from the evaluation of the free 
energy of a trial state of the form (13.91) . Using some projected states in an appropriate way, 
this can be done by adapting well-known tools of finite dimensional semi-classical analysis. 
The lower bound is as usual the tougher part, and we will try to give (some of) the main 
ideas in the sequel. 



F c i[o] = J F NL [u\do(u) + J ^-log 


4. Some tools of the proof 

The first point is to characterize the (limits of) quantum states over the Fock space 
in terms of classical objects, i.e. probability measures over the one-body Hilbert space. 
This is the object of the so-called quantum de Finetti theorem, that we discuss first. It 
will be rather clear from this kind of results how one can deal with the interaction energy 
term in (13.111) . The relative entropy term requires new specific tools, discussed in a second 
subsection. 

4.1. Quantum de Finetti measures. The main message of the quantum de Finetti the¬ 
orem is informally that, “For many practical purposes, any bosonic state looks like a convex 
combination of product states”. A little bit more precisely, if one picks a bosonic A-body 
state, with A large, and look at its A:-body density matrix, for k small, one obtains “almost” 
a convex combination of product states. 

There is a vast literature devoted to rigorous statements and proofs of these claims. We 
shall not review it here in details, see [mm si- For the purpose of this note, we rely 
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on recent approaches and versions of the theorem due mainly to Ammari-Nier and Lewin- 
Nam-Rougerie 12113, HI 128]. To deal with the case where p > 1 in (12.711 . some improvements 
are necessary, provided in m, but for the proof of the previous theorems, the following 
result of [2] is sufficient 


Theorem 4.1 (Grand-canonical quantum de Finetti theorem). 

Let 0 ^ T„ be a sequence of states on the Fock space Assume that there exists a 

sequence 0 < e n —>■ 0 such that, for all k £ N 


etTr, 


p(fc) 


^ Ck < oo. 


(4.1) 


Then, there exists a unique Borel probability measure u on S) such that, along a subsequence, 

k\{e n ) k T^ ^ [ \u® k )(u® k \dv(u) (4.2) 

■’s) 

weakly in the trace-class for every integer k. 


Some comments: 

• The measure v is called the de Finetti measure of the sequence T n . It is the classical 
object we are after. Its existence is a rather general fact, based mostly on the bosonic 
symmetry of the states under consideration. Theorems 13.II and 13.21 can be rephrased 
as saying that the de Finetti measure of the sequence Tyr is the nonlinear Gibbs 
measure p. 

• Here e n should be thought of as a semiclassical parameter essentially forcing the 
particle number to behave as e” 1 . To apply this result in our context, we need to 
check that m holds, with e n = T 1 and T = Ty^. In the context of this note, 
this is not too difficult because the particle number is well under control. Bounds 
on the expectation of Af k in some state naturally translate into trace-class bounds 
on the k -body density matrix. 

• Since we here get convergence of density matrices, a lower bound to the interaction 
energy term in (13.111) follows from this theorem. 


4.2. Berezin-Lieb inequality for the relative entropy. There remains to understand 
how to estimate the relative entropy term in (13.111) . One should keep the following in mind: 

• Contrarily to the energy, it is not expressed in terms of reduced density matrices 
only. It genuinely depends of the full state. 

• One should not “undo” the relative entropy to consider only the von Neunmann 
entropies of To,t and T\ t directly, since we do not expect a good control on those. 

The main tool to link quantum entropy terms to their semiclassical counterparts are the 
so-called Berezin-Lieb inequalities EES US]. Since we work in infinite dimensions and we 
cannot undo the relative entropy to rely on known inequalities, we need the following new 
result: 


Theorem 4.2 (Berezin-Lieb inequality for the relative entropy). 

Let e n —>• 0 and {T n } ; {T^} be two sequences of states satisfying the assumptions of Theo¬ 
rem \4-1\ Let u and v' be their de Finetti measures. Then 

“SS 1 ^ A ( logr “ - logI t )] 58 /f? log (f?) dl/ ' 


(4.3) 
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Some comments: 

• Clearly, this will give the desired lower bound to the relative entropy term in (13.111) 
if we can prove first that the de Finetti measure of the free Gibbs state Fo,t is given 
by hq. Since all the density matrices of Fo,t can be explicitly computed, this is 
easily shown by direct inspection. 

• For the proof of this result, we rely heavily on a constructive approach to the proof 
of Theorem rm Indeed, not only the existence of the measure is crucial, but also 
the precise way one can approach it by explicit functions of the states. 

• The link with semi-classics proceeds by seeing the de Finetti measure both as a 
lower symbol/covariant symbol/Husimi function/anti-Wick quantization and as an 
approximate upper symbol/contravariant symbol/Wigner measure/Wick quantiza¬ 
tion based on a coherent state decomposition (made of projectors onto states of the 
form (13.101) 1. 

• We also rely on deep properties of the quantum relative entropy: joint convexity 
and monotonicity under two-positive trace preserving maps mmm- These are 
in fact equivalent to the strong subadditivity of the quantum entropy, a crucial 
property proved by Lieb and Ruskai [Ml El]. 
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